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O ■ Abstract 

> 

^ , Matrix Riccati equations and other nonlinear ordinary differential equations 

H : 

>y , with superposition formulas are, in the case of constant coefficients, shown 

to have the same exact solutions as their group theoretical discretizations. 
Explicit solutions of certain classes of scalar and matrix Riccati equations are 
presented as an illustration of the general results. 
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Typeset using REVT^ 



I. INTRODUCTION 



The concept of linear superposition of solutions of a linear equation, the Schrodinger 
equation, forms the basis of all of quantum mechanics. It has two different aspects. One is 
that a linear combination of two solutions is again a solution. The second is that the general 
solution of a linear equation can be written as a linear combination of a sufficient number 
of particular solutions. 

Both of these aspects of superposition have their counterpart in the theory of nonlinear 
phenomena, as represented by nonlinear differential and difference equations. Thus, in soli- 
ton theory, specific solutions, namely solitons, can be "superposed" in a nonlinear manner, 
to obtain multisoliton solutions of integrable nonlinear partial differential equations . On 
the other hand, certain classes of nonlinear ordinary equations exist, that allow a nonlinear 
superposition formula, in the sense that their general solution can be expressed analytically 
in terms of an appropriate number of particular solutions. 

These two different aspects of "nonlinear superposition" are actually related via Backlund 
t r ansf ormat ions . 

Several recent articles have been devoted to a study of nonlinear ordinary differential 
equations (ODEs) with superposition formulas 0-0- The entire approach is based on a 
theorem due to Sophus Lie that states the following. 

Theorem 1 (S. Lie) The general solution y{t; Ci, . . . , c„) of a system of n ordinary differ- 
ential equations 

V=v(y,t), yeC'^forR^) (i.i) 

can be written as a function of m particular solutions and n arbitrary constants 

y = F{yi{t),...,ymit);ci,...,Cn)) (1.2) 

if and only if: 
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1. The system has the form 



V{t) = J2 Zk{t)l]k{y, t), 1 < r < oo (1.3) 

A;=l 

2. The vector fields 

n 

Xk = Y.^mdv.^ (1-4) 

generate a Lie algebra L of finite dimension dimL = r < oo. The number of equations 
n, number of fundamental solutions m and the dimension of the Lie algebra r satisfy 

nm > r. (1-5) 

A geometric approach to such equations has been developed p. The variables yi,. ■ ■ , yn 
are considered to be coordinates on a homogeneous space M ~ G/ Gq with the group G acting 
transitively on M. The Lie algebra L of vector fields ( |1.4| ) is the Lie algebra corresponding 
to the group G. The subgroup Gq is the isotropy group of the origin and its Lie algebra 
Lq is the subalgebra of L consisting of vector fields that vanish at the origin. Thus, it is 
possible to associate a system of ODEs satisfying Lie's criterion with every group-subgroup 
pair (Go, G), or with every Lie algebra pair (Lq, L), Lq G L. 

We call the system (|1.3|) decomposable if a subsystem of equations can be split off that 
involves only a subset of variables, namely, yi,. . . , yk, I < k < n, that itself satisfies Lie's 
theorem. This occurs if the action of G on M is not primitive, i.e., if a G-invariant foliation 
of M exists Indecomposable systems of ODE's with superposition formulas are obtained if 
the pair of Lie algebras (Lq, L) forms a transitive primitive Lie algebra P-p!3[]. This requires 



that Lq be a maximal subalgebra of L, not containing any ideal of L. Transitive primitive 
Lie algebras over C and M have been classified P-p!3[]. Let us now consider a particularly 
interesting class of transitive primitive Lie algebras, namely one where L is a simple classical 
Lie algebra over C or R and Lq is a maximal parabolic subalgebra. The space M ~ G/Gq 
is then a Grassmanian. Using affine coordinates on this space we obtain vector fields with 



polynomial coefficients of order at most four P|. The corresponding ODEs (|1.3| ) have the 
same polynomial nonlinearities. 

These equations have the following properties: 

1. By construction, they allow a nonlinear superposition formula. 

2. They are linearizable by introducing homogeneous coordinates on the Grassmannian, 
i.e., by an embedding into a larger system of linear equations. 

3. They have the Painleve property [||], i.e., their singularity structure is particularly 
simple: the only movable singularities of their solutions are poles. 

A discretization of the nonlinear ODEs with superposition formulas has been proposed in 



a recent article |14|. It preserves all of these properties, in particular the Painleve property, 



manifesting itself as singularity confinement ||T5 



The discretization was performed for the case when L was any classical complex Lie 
algebra and Lq any of its maximal parabolic subalgebras. The discretization made use of 
the linearization of the equations ( |1.3| ). It can, however, be presented schematically without 
using the corresponding linear system. 

Indeed, let ( [1.4|) represent the vector fields of Lie's theorem. Eq. (|1.3| ) can be rewritten 



as 



y = X{t)y (1.6) 

r n 

x{t) = j2z,{t)j2vtm., (1-7) 

k=i fj.=i 

As t varies, the vector field X(t) follows some path in the Lie algebra L. The vector 
fields ( p..4| ) can be integrated to give the corresponding local group action 

y = gu, (1.8) 

where u is some chosen point on M. The form of action ( |1.8| ) is independent of the equation 
( |1.6|) , i.e. of the functions Zk{t). It involves r = dimL = dimC group parameters. The 



general solution of eq. ( |1.6|) can be represented in the form (|1.8|) , where g = g{t) is a function 
of t. As t varies, g{t) will trace a path in the group manifold, corresponding to the path 
that X{t) follows in the Lie algebra. Reconstructing this path is equivalent to solving the 
equation (|1.6|). 

The difference equations with all the desired properties are obtained in the form ( |1.8[ ) 
with g a function of time (i.e., each of the r group parameters in g depends on time). The 
equations are 

y = 9y, y = y{t), y = y{t + a), (1.9) 

where a is the spacing of the lattice. The group element g = g(t) is to be evaluated at the 
same point as y{t) on the right hand side. For a — we have g{t) = I + aX(t) and the 
difference equation ( p.. 91 ) goes into the differential equation ( p.. 6] ). 

The purpose of this article is to explore some properties of the difference equations ( |1.9[ ) 
and in particular to show that in the case of constant coefficients in the differential equations 
( |L6|) the solutions of the continuous and discrete equations coincide exactly. 

A general theorem on solutions of continuous and discrete matrix Riccati equation is 
formulated and proven in Section O. Section 111 is devoted to examples. 



II. CONTINUOUS AND DISCRETE MATRIX RICCATI EQUATIONS 



Let us consider the Grassmannian M ~ G/Gq where G is SL(A^, C) and Go is a maximal 
parabolic subgroup. We have 



N = p + q, G~ 



Lril Ct12 



G2I G22 



Gq 



^ Gn ^ 



G21 G22 



(2.1) 



where Gn G C^^^ G12 G Cp^'?,G2i G C«^^ G22 G O^" and M = pg, detG = 1. Affine 
coordinates on this space (in the neighborhood of the origin) are given by the matrix elements 
of a matrix W G C^^''. The nonlinear ODEs with superposition formulas in this case are 
the Matrix Riccati Equations ||^-[7|,p!6[ 



W = A + BW + WC + WDW, 



(2.2) 



where 



^ B{t) A{t) ^ 



V 



, TiB = TrC, 



(2.3) 



-m -c{t) 

For t fixed, ^(t) is an element of tlie Lie algebra sl(A^, C). Tlie division into blocks is tlie 
same in eq. (|2.3|) as in (p.l|) Tlie matrices A, B, C, D are known functions of t. As t varies 
tlie matrix ^{t) follows some path lying in sl(A^, C). 

The Matrix Riccati Equation ( p.2|) is linearized by introducing homogeneous coordinates, 
i.e. the matrix elements of two matrices X and Y . We put 



W = XY-\ X e C^^^ Y e a'"'. 



(2.4) 



The origin on M corresponds to X = O, Y = I. When the matrices X, Y satisfy the linear 
system 



\ 



Bit) A(t) 

-Bit) -at) 



the matrix W = XY ^ satisfies eq. (|2^) . 

The ordinary Riccati equation for w{t) a scalar function, namely 

w = a + {b + c)w + dw'^, 



(2.5) 



(2.6) 



is obtained for p = g = 1. For p > 2, g = 1 eq. (|2.2| ) is called the projective Riccati equation, 
since it corresponds to the projective action of SL(p + 1, C) on the space C^. 
The discretization of eq. (|2.5|) is 



(2.7) 



Bit) A(t) 

y -Bit) -c(t) J 

where X = X{t + a), Y = Y{t + a) and the discrete matrix Riccati equation is obtained via 
the projection (|2.4| ), i.e.. 



W = {GiiW + Gi2){G2iW + G22)~^ . (2.8) 

Now let us restrict ourselves to the case of the Matrix Riccati Equations ( |2.2| ) with con- 
stant coefficients. In this case the matrix ^ of eq. ( p.3| ) can be viewed as a one-dimensional 
subalgebra of sl(A^, C) and the corresponding one-dimensional subgroup of SL(A^, C) is ob- 
tained by exponentiating ^ 

G = e< = I + ei + ^e^e + ---- (2.9) 



This matrix G is to be inserted into eq. (|2.71 ) and eq. ( p.8|) , first identifying the group 
parameter e with the lattice spacing a, for example, putting e = a. Alternatively the form 
of the transformation G can be obtained by integrating the vector fields on the right hand 
side of eq. ( p.2|) . To do this we rewrite eq. ( |2.2|) in components as 

W^a = A^a + Bo^pWpa + W^a + W^bD.pWpa, (2.10) 

where l<a<p, l<a<q' (repeated Greek labels are summed from 1 to p, Latin ones 
from 1 to g). Eq. (|2.8|) is obtained directly by integrating 

= A^a + B^pWpa + W^,G,a + W ^l^D^pW (2.11) 

where Wc,a\e=0 = W^a- 

We now come to our main result which we formulate as a theorem. 



Theorem 2 The general solution of the continuous matrix Riccati equation ( \2.3i ) with con- 
stant coefficients coincides with the general solution of the discrete matrix Riccati equation 
^7^), with coefficients G^^ satisfying eq. ([17 



Proof. Let us rewrite the discrete linear system (|2.7| ) as 

(2.12) 



and expand both sides in powers of a: 



Comparing coefficients of the powers of a we have 

Jn 

—Z = CZ. (2.13) 



For n = this is satisfied trivially, for n = 1 we obtain the continuous linear system 

Z = iZ. (2.14) 



Let Z{t) be a solution of eq. ( 2.14 ) and hence W a solution of the continuous matrix Riccati 
equation (|2.2| ). In order for Z to be a solution of the discrete system ( p.7|) , and hence W a 
solution of the discrete matrix Riccati equation (|2.8|) , Z must also satisfy eq. (|2.13|) for all 
77, > 2. For n = 2 we have 

^ = #(e^) = ?^ + e^ = e^^> (2-15) 
dt^ dV^ ' dt dt ^ ' ^ ' 

since d^jdt = ^ and Z = ^Z. For > 3 eq. (|2.13|) follows by induction, hence every solution 
of eq. (|2.14| ) with constant ^ satisfies the difference equation ([2171 ). □ 



The general solution of eq. ( p.5|) depends on p.q significant free parameters, the initial 
conditions. The same is true for the difference equation (|2.?1 ). Hence, the general solution of 
the differential equation ( |2.5|) provides the general solution of the difference equation ( p.7|) 
and vice versa. 

III. EXAMPLES 

A. The Scalar Riccati Equation 

Let us ffist consider a scalar Riccati equation with real constant coefficients: 

u = a + bu + cu'^, c 7^ 0, (3.1) 

and introduce A = 6^ — 4ac. We simplify this equation by using an appropriate linear 
transformation u = ay + (3. Depending on the sign of A we obtain three different standard 
forms of (|3.1|) . 



1. A>0 



We choose a and /3 so as to obtain 

y = uj{y'^ — 1), ^ = const. (3.2) 

Integrating the vector field X = (y^ + l)d/dy we obtain the discrete Riccati equation 

y cosh auj — sinh au . , 

y = ^ ' ^-3) 

—y smh auj + cosh auj 

where y = y{t), y = y(t + a). It is easy to verify by direct calculation that the general 
solution of both equations (|3^ ) , (|3^ ) is 



y = - cothLj{t-to), (3.4) 
where to is an arbitrary constant. 

2. A<0 

An appropriate choice of a and P yields 

y = ojiy^ + 1), ^ = const, (3.5) 

Integration of the vector field X = {y'^ — l)d/dy provides us with the discrete Riccati equation 

y cos au + sin au . . 

y = : . (3.6) 

—y sm au + cos au 

The general solution of these equations is 

y = tan Lj{t — to), (3.7) 

where to is an arbitrary constant. 
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5. A = 



In this case the Riccati equation (|3.1|) can be reduced to 

y = y\ (3.8) 

and its discrete form is 

y = -y—. (3.9) 

1 - ay 



The general solution of both equations ( p.8|) and (|3.9| ) is 



y = -T^^ (3-10) 

t — to 



where to is an arbitrary constant. 



B. Coupled Riccati Equations Based on a Nonprimitive Group Action 

Most of the earlier articles on nonlinear ODEs with superposition formulas concerned in- 
decomposable systems of equations. A study of imprimitive actions and hence decomposable 
systems was initiated quite recently |T7| . 

Here we shall give an example of such a system and its discretization. Consider the Lie 
algebra sl(2,M), realized as a subalgebra of diff(2,M), by the vector fields. 

Xi = <9,, X2 = xd^+ydy, Xs = x^d, + {2xy + ky'')dy, (3.11) 

when k = 0, or k = 1. 

The corresponding ODEs are 

x = a + bx + cx'^, y = by + c{2xy + ky'^) , (3.12) 

Let us now restrict to the case when a, b, c are constants and moreover assume that b'^—Aac < 
0, c 7^ 0. We can then, by means of a linear transformation of x take eq. ( ^.12| ) into the 
standard form 
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x = uj{x'^ + l), y = uj{2xy + ky'^), = const. (3.13) 

To obtain a discretization that preserves the underlying group structure, we integrate the 
vector field 

X = + x^)d, + {2xy + ky^)dyl (3.14) 



We have 



'^ = u{l + x% ^=uj{2xy + ky% (3.15) 



where x\u=o = x, y\a=o = y- The result is 

X cos au + sin au 



X = : , 



yix"^ + 1) , , 

y= ^^ ^ \ (3.16) 



—X sm auj + cos au 
D 

where D = {x^ + 1 + kxy) (— x sin au + cos auY — ky{x cos au + sin au) {—x sin au + cos au) 



Thus, eq. ( p. 161 ) is a difference equation that has eq. ( p.l3| ) as its continuous limit (for 
a 0). 



Eq. (|3.13| ) is easy to solve and we obtain 



X = tanti;(t — to), 

^ a;cos^ci;(t — to) ~ ^[sin2ti;(t — to)]/2' ^ ^ 



where a and to are integration constants. Expression ( p.l7|) is also the general solution of 
the difference equation ( p.l6|) . 



C. Projective Riccati Equations 



Let us now consider the matrix Riccati equation ( |2.2| ) for q = 1 which is called the 
projective Riccati equations. 
We rewrite the system as 
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y = a + By + y^Cy 



(3.18) 



where a, y G are column vectors, B, C G /^pxp matrices and T denotes transposition. 
We consider the case when all coefficients are constants. 

We simplify eq. (|3.18|) by a linear transformation with constant coefficients 



y = Au + p (3.19) 

Choosing A and p appropriately we can, for instance, transform C and B in such a way that 

a = 5,i, % = 0, j = l,2,...,p (3.20) 

with some further freedom remaining. 

Now let us consider the simplest case of p = 2. Changing notations we write the trans- 
formed projective Riccati equation as 

X = 1 + y + x'^, y = a + by + xy, (3-21) 

where a, b are constants. 

In order to discretize the equations we must either exponentiate a matrix as in eq. ( p.9|) 
or integrate the vector field 

X = {l + y + x'')d^ + {a + by + xy)dy, (3.22) 

We choose the second procedure and put 

^ = {l+y + x'), ^ = {a + by + xy), (3.23) 
da da 

where x\a=o = x, y\a=o = y. 

Eliminating y from eq. (|3.23|) we have 

X = 3xx + bx — x"^ — bxP' — X + a — b (3.24) 

where dots denote differentiation with respect to a. 
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Eq. (p.24|) has the Painleve property and can be transformed to one of the 50 stan- 



dard forms 1 18 1, studied by Painleve and Gambier. More specifically it is equivalent to the 



linearizable equation No. VI in the list reproduced by Ince [|18| (see p. 334). Indeed, we put 

x{a) = a{cx)w{T) + l3, r = t{cx), (3.25) 

while the constant (3 must satisfy: 

(3^ + b(3'^ + (3 + b- a = 

The functions a{a) and T{a) must satisfy linear equations with constant coefficients: 

a-{b + 3p)a + (3/?^ + 2bp + l)a = (3.26) 
f = -a. (3.27) 

The function w(t) then satisfies 



Wt-t = —•^wWr — w + q{T){wT- + w ) (3.28) 

3d — 3a3 — ba , , 

q{r) = . (3.29) 



Equation (|3.28|) is linearized by the Cole-Hopf transformation 



fr 

W = y, frrr = q{'r)frr- (3.30) 

Eq. (|3.30|) can be solved in quadratures, however, we shall only consider the simplest case, 
when g(r) = . This implies a constraint between the constants a, 6, namely 

,M.O, (3,31) 



Equation (|3.30| ) in this case is solved trivially. Returning to the original variables x and y 
we have 

^ 2Pe^- + g 
^ ' Pe^" + Q + Re-'"' 

^("^ = ^ Pe^^ + Q + i?e-^^ ^ ^^-^2^ 
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Here P, Q, R are integration constants and we have put 



T = e 



k(7 



a = —ke 



k(T 



(3.33) 



We ehminate the constants P,Q,R by imposing initial conditions as in eq. ( |3.23| ). Finally, 
we obtain difference equations 



where 



X 



9iix + guy + 913 

gsix + g32y + gss ' 



y 



g2ix + g22y + g23 
gsix + g32y + gss ' 



(3.34) 



g^^ = {p- 2k){2p - A;)e^" - 4(/3 - kf + (3{2(3 - 3k)e~'"', 

gi2 = {P- 2k)e''^ + 2{k -(3)+ pe-'"', 

^13 = {(3- 2k){-l3^ + kp + l)e'^'^ - 2{p -k)+ pe-'"' 

g2i = {Akp - Ak^ -/3^- l){2/3 - A;)e^" - A{2f3k - k^ - (3^ - 1){(3 - k) 
-(/32 + l)(2/5-3A;)e-'=^ 
= {Akp - -f3^- l)e'=" - 2{2pk - k^ - - 1) - {(3^ + 1)6"'=", 
= [Akf] - Ae -f3^- l){-f3^ + f3k + l)e^" - 2{2(3k - ^ - f3^ - l){2(3k + 1) 
+ {(3^ + l)(/?2 + 2e - 3(3k + l)e~'^", 

^31 = (2/5 - k)e''' - 4(/? -k) + (2/5 - 3A;)e-^^ 

g,, = e'^-2 + e-'% 



gs3 = {l+(3k- (3')e'''' + 2{(3' - 2(3k + 1) + {-(3' - 2k' + ?,(3k + l)e 



'k(j 



(3.35) 



As in the previous examples the general solutions of the projective Riccati equations ( p. 211) 
and their discretization ( p.34| ), coincide. 



IV. CONCLUSION 

The results of this article can be summed up as follows. 
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1. We have shown that matrix Riccati equations with constant coefficients and their 
group theoretical discretizations have the same solutions. 

2. We have shown how other types of equations with superposition formulas can be 
discretized via an integration of the underlying vector fields. For constant coefficients 
the continuous and discrete equations again have the same solutions. 

An open question, presently under study, is whether the results of this article can be used 
to obtain an effective perturbative approach to solving differential and difference equations 
that allow superposition formulas, but do not have constant coefficients. 
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